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Abstract 

Within the four-particle cluster approximation for the proton 
ordering model we study the influence of shear stress eg on the 
phase transition, static dielectric, elastic and thermal properties of 
deuterated KD2PO4 type ferroelectrics. Thermodynamic poten- 
tials and physical characteristics of the crystals in the presence 
of stress as are calculated. Numerical analysis of the obtained 
results is performed, and the set of the theory parameters provid- 
ing the best fit to the available experimental data is found. The 
To — as phase diagram is constructed; temperature dependences of 
the calculated quantities are studied, possible stress dependences 
are indicated. Influence of as stress is analogous to that of electric 
field conjugate to the order parameter. 
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1 Introduction 



More that sixty years have passed since the discovery of the ferroelec- 
tricity in the KH2PO4 crystals. During those years, a huge number of 
papers has been devoted to the studies of the phase transitions in the 
crystals of the KH2PO4 family and of their physical properties. The 
most prominent peculiarity of these studies is a strong connection be- 
tween theory and experiment, which is thought to be an important source 
of the obtained progress in microscopic understanding of their proper- 
ties. In this aspect, the high pressure studies of these crystals come 
extremely important. External pressure changes the internal structure 
of the system, altering thereby the molecular potentials in the crystal 
and its dielectric and thermal properties. And shear stresses give rise to 
shear piezoelectric lattice strains induced also by external electric fields 
via the piezoelectric effect; that allows one to consistently explore the 
electric, electromechanic and thermal characteristics of the crystals. 

A microscopic model of strained KH2PO4 type crystals has been pro- 
posed in jy, g, ||. According to this model, external mechanical stresses 
give rise to additional internal fields, linear in strains and (in the case of 
diagonal components of stress tensor) mean values of quasispins. Influ- 
ence of stresses of different symmetries on energies of dcutcron configu- 
rations was studied. 

In § || there were also explored the effects of ayi = &xx — &yy 
stress on the transition temperature to the ferroelectric phase and the 
relations between the applied stress a xx — a yy and the induced strain 
£12 = £xx — £yy taking into account the deuteron rearrangement on hy- 
drogen bonds. In g § using the model proposed in jy, [|, within 
the cluster approximation taking into account the short-range and long- 
range interactions and stress (T12, the dielectric, piezoelectric, and ther- 
mal characteristics of KD2PO4 were calculated and studied. It has been 
shown that the stress tri2 can lead to a phase transition to the monoclinic 
phase. 

Influence of hydrostatic pressure and uniaxial pressure — p — 03 on 
the physical properties of the KD2PO4 type crystals was investigated 
in H 0, [§. We showed that under the proper choice of the theory 
parameters, a satisfactory description of the experimental data for the 
pressure dependences of spontaneous polarization, longitudinal static di- 
electric permittivity and transition temperature was obtained. Using the 
Glauber model j|] the expressions for the real and imaginary parts of 
the longitudinal dynamic dielectric permittivity of KD2PO4 were found, 
and their temperature and frequency dependences at different values of 
hydrostatic pressure were calculated. 

It is also interesting to study the influence of the shear stress 0% giving 
rise to a shear strain £6 — £xy, which transforms after the irreducible 



1 



representation Bi. After this representation, polarization P3 transforms 
too; hence, the influence of external stress o§ = a xy is similar to that of 
electric field E3. Importance of this study results also from the fact that 
spontaneous polarization P3 in KD2PO4 is accompanied by spontaneous 
shear strain eg. 

In [jlO| the dielectric, piezoelectric, and elastic properties of KH2PO4 
with taking into account the strain Eq induced by the piezoelectric effect 
are studied within the modified Slater model. 

Experimental measurements of dielectric, electromechanic, and elas- 
tic characteristics of K^Hi-^D^PC^ were performed in several works. 
Thus, in the temperature dependences of the dielectric permittivities 
of free and clamped KH2PO4 crystal, piezoelectric constants d^e, e36, /136 
and elastic constants Cq 6 , c§ 6 and Sq 6 , Sq 6 were reported. Temperature 
dependences of the piezoelectric module d^a of the KH2PO4 crystal at 
direct and inverse piezoelectric effect are given in jl^] and JD|, respec- 
tively. Temperature dependence of elastic constant cf 6 for KH2PO4 is 
presented in [jl4|. Experimental data for the temperature dependence 
of g?36, £33 and Sq & for a partially deuterated crystal with x = 0.89 are 
given in |l5|, 

In this paper we study the influence of the shear stress <j§ on the phase 
transition and physical properties of deuterated ferroelectric crystals of 
the KH2PO4 family. We shall also compare the obtained results for 
thermodynamic, dielectric, elastic, and piezoelectric characteristics of 
the crystals with the corresponding experimental data. 



2 The crystal Hamiltonian 

We consider a system of deuterons moving on -D-. . .- bonds in deuter- 
ated crystals of the KD2PO4 type. The primitive cell of such a crystal 
is composed of two neighbouring PO4 tetrahedra together with four hy- 
drogen bonds attached to one of them ("A" type tetrahedra). Hydrogen 
bonds going to another ("£?" type) tetrahedron belong to four nearest 
structural elements surrounding it (see the figure below). 

Hamiltonian of a deuteron subsystem of the KD2PO4 type crystals 
taking into account the short-range and long-range interactions in the 
presence of the shear stress o§ = <j xy giving rise to the strain e§ — 
e xy and of the external electric fields Ei{i — 1,2,3) applied along the 
crystallographic axes (a, 6, c) reads 

Hi = ^cgfeg - vNe° 36 E 3 e e - + (2.1) 

q'f'qf qf 
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Figure 1: Primitive cell of the KD2PO4 crystal. Numbers in circles correspond 
to hydrogen bonds; 1,2 are the deuteron equilibrium sites. By dashed lines 
we show the directions of hydrogen bonds in an unstrained crystal in the 
paraelectric phase, a, b are the axes of an undeformed 742d cell; 7 is an angle 
between the directions of hydrogen bonds in strained and unstrained crystals. 



The first term in the right hand side of (2.1) corresponds to that part 
of the elastic energy which does not depend on hydrogen arrangement 
(c§q is the so called "seed" elastic constant); the second term in (|2.l[ ) 
is the energy of interaction between polarization induced by the piezo- 
electric effect due to the strain e 6 (not taking into account the hydrogen 
subsystem contribution) and field £3 (e§ 6 is the seed coefficient of the 
piezoelectric stress); the third term is the energy due to polarization in- 
duced by external field independently of the hydrogen configuration (^33 



is the seed dielectric susceptibility). The last term in (2.1) describes the 
short range configurational interactions between deuterons around the 
"A" and "B" type tetrahedra. Two eigenvalues of Ising spin a q f = ±1 
are assigned to two equilibrium positions of a deuteron on the /-th bond 
in the q-th unit cell (a q t — ±1). v — f/fce, v is the primitive cell vol- 
ume; fee is the Boltzmann constant. Fl are internal fields, arising in 
strained crystals, which include the effective long-range interaction be- 
tween deuterons (taken into account in the mean field approximation) 
and the additional internal fields related to the strain Sq ||: 

2fiF z (6) = 2/xi^(6) = -2/zF|(6) = -2/iF|(6) = 2/xF|(6) = 

= 2i/ c 7/ 1)z (6) - 2V> 6 £6, 
2/^(6) = 2v lV [ 1)x - v (6) + 2u 3 4 1)x ' v (6) + 

+ 2z, 2 [^(6)+7 7 i 1 ^(6)]-2^ 6£6 , 
2^(6) = 2 I rf^(6) + ^*( 6 )] + 
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+ 2v 1 4 1)x ' v (6) + 2v 3 r,P x ' v (6) - 20 6 £ 6 , (2.2) 

"(6) = 2v 3V [ 1)x > y (6) + 2v 1 r 1 { 3 1)x ' y (6) + 
+ 2Mv { 2 1)x ' y (S) + V { l )x ' y m-^e 6 , 
2^(6) = 2 I rf^(6) + 4 1)a: ' y (6)] + 
+ 2v 3V i 1 ^(6) + 2v lV [ 1 ^(6)-2^e 6 , 

where 

^(6) = <a gl )* = (a q2 y = (a q3 y = <<7, 4 )*, ^ 1 >-»(6) = (cr,/)-*; 

. „ Jll ^12 ^13 

^ c = z/i + 2i/ 2 + v z \ v\ = —,v 2 = —,v 3 = — ; 

^//' = ZIr -r , Jff'ilQ') ls the Fourier transform of the long-range 
deuteron-deuteron interaction, -06 is the so called deformation potential; 
H = eS is the dipole moment of a hydrogen bond; 8 is the D-site distance. 



In (2.1) and (2.2) we do not take into account piezoelectric shear 
strains £4 or £5 induced by electric fields E± and E2, respectively. 

The Hamiltonian is usually chosen such that to reproduce the energy 
levels of the Slater- type model for KDP (see, for instance fl7| ) - the 
Slater energies £, w, and aii (e < i» C Wx), determined by the energies of 
up-down £ s , lateral e a , single-ionized e%, and double-ionized £0 deuteron 
configurations. 

£ = £ a -£ s , w = si—e s , wi=e -e s . 

If £6 = 0, configurations with two hydrogens in potential wells being 
close to upper and lower oxygens of a given PO4 group (each PO4 tetra- 
hedron is oriented such that two of its edges are parallel to the ab plane) 
and with the hydrogens on the two other bonds being close to the neigh- 
boring tetrahedra, have the same energy e s , assumed to be the lowest. 
Correspondingly, lateral configurations with two hydrogens close to an 
upper and a lower oxygens are four-fold degenerated; single-ionized with 
only one (or three) hydrogens close to a given group are eight-fold de- 
generated, and double-ionized with four hydrogens (or without any) are 
twice degenerated. 

However, the strain <j§ splits certain deuteron configurations of the 
conventional Slater- Takagi model. In Table below we present all pos- 
sible deuteron configurations and their energies. Since the system is 
not any longer symmetric with respect to a reflection in the ab plane 
<Th or to the reflection with 7r/4 rotation around the c-axis S4 (both 
operations change the sign of polarization and strain £g), the up-down 
configurations (i = 1,2) split to two different levels, and lateral config- 
urations split to two groups of twice degenerated levels each (i = 5,6) 
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and (i = 7,8). Configurations within each group are symmetric with 
respect to it/2 rotation C%. Similarly, single- ionized configuration are 
divided into two groups, within each the c-component of dipole mo- 
ments assigned to a configuration are directed up (i = 9, 10, 11, 12) or 
down (i = 13,14,15,16). 

Here we assume that the strain Eq alters the energies of deuteron 
configurations only by splitting the degenerated levels due to lowering 
of the system symmetry. From geometric point of view, this happens 
mostly because the angle between perpendicular in the paraelectric phase 
hydrogen bonds is changed, and the PO4 groups are distorted, whereas 
in the hydrostatic pressure case the changes are usually attributed to 
pressure-induced changes in the D-site distance S || [ll| . 

To rewrite the energies of deuteron configurations in terms of 
pseudospins, we associate the configuration operator Ni with the con- 
figuration i according to the following rule: each operator is a product 
of four factors, one per each hydrogen bond, each factor being equal to 
o(l + c<7/) if deuteron is in the first minimum at the f-th bond and 
1(1 — a qf) otherwise. Then the Hamiltonian of the short-range interac- 
tions takes the form 

N 16 

4 hort (6) =Y,Y,[Nf{q) +Nf (q)] = (2.3) 

q=l i=l 
JV f 4 

9=1 I /=1 



Vql Vq2 & q3 & ql & ql & q& ' q\ ° q3 °\}4 & q2 & qi O ' ? 4 

~2 2 2~ ~2 2 2~ ~2 2 2~ ~2 2 2~ 



+ (V + 6 6 e 6 ) "-f + +(V- 6 6 e e ) + £«L 



2 2 2 2 



U 



C ql & q3 . & q2 & qi 

~2 2~ ~2 2~ 



L 2 2 2 2 

N 16 



q=l i=l 



Here the following notations are used 

V = V13 = V 2 3 = V34 = V41 = 

U = Vi3 = V 2 a = tjWi - £, <f> = 4e + 2wi - 8w, 

E = E a -E s , W = ei-E s , Wi=£o-£ s , 

where e 8 , e a , £1, £0 are energies of up-down, lateral, single-ionized and 
double-ionized deuteron configurations in an unstrained crystal, respec- 
tively. 
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The further calculations will be performedd in the four-particle clus- 
ter approximation. The cluster Hamiltonian in the presence of the stress 
erg and electric field Ei reads 



^4(6) 



(s s6 + 2s 16 ) y: a -f 



(2.4) 



-£6(5s6+25ie) 
f {V + S 6 e 6 ) 



Cgl &q2 C<j3 Cgl °92 Cg4 

~2 2 2~ ~2 2 2~ 

Cql <J q2 Cq3 &qi 

~2 2~ ~2 2~ 



fg3 &q4 CT g 2 0" g 3 CT g 4 

~2 2 2~ ~2 2 2~ 
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~2 2~ 



Cq2 0" g 4 

~2 2~ 



^Ogi 0^2 

2 2 2 2 



E 



Z /6 gg/ 

/? 2 ' 



In (2.4) 



6 = *?6 = *26 = 4i - 4» = /3[- A 6 + 2^(6) + Ma^]; 

f 6 = fl-Af + 2/*f?(6) + Mx cos 7 i?i], 
f 6 = /J[-Ag + 2^(6) - Mx sin 7^], 
/3[-A| + 2 M if(6)- M ± cos 7 i?i], 

/3[-Ag + 2Aif?(6)+At±sin7£7i]; (2.5) 
/?[-Ag + 2/^(6) + Mx sin 7j B 2 ], 
/?[-Ag + 2^1(6) - /x_L cos 7 S 2 ], 
/3[-A 6 y + 2 M F|(6) - M x sin 7j B 2 ], 
/?[-A* + 2 M F|(6) + /x_L cos 7^2]; 



'36 



'46 



-1/ _ 

'36 — 



^46 

/Lt3 and /Ltj_ are the longitudinal and transverse effective dipole moments 
of primitive cells, created by displacements and polarization of heavy 
ions which are triggered by deuteron ordering. 

We took into account the following relations which are obeyed in the 
presence of stress and fields Ei 



Ml3 = A*23 = M33 = M43 = M3; 

Mil = Mi cos 7, /! 2 i = -fix sin 7, 
M3i = -MJ-COS7, fin = fi± sin 7 

Using the condition 



M12 = sin 7, 
M32 = -MX sin 7, 



(2.6) 

-M22 = MX cos 7; 
M42 = Mx cos 7. 



s P K /e -^ 6 >} 

Sp{e-^< 6 )} S P {e-^/( 6 )} ' 

where the single-particle deuteron Hamiltonians are 



fyi ( 6 ) = 



4 /6 



*/6> 



(2.7) 



(2.8) 



7 



we calculate the single-particle distribution functions of deuterons and 
exclude the parameters A 6 



^(6) = ^, 



(2.9) 



where 



m z {6) = sinh(2^ + (3S s6 s 6 ) + 26sinh(^ - (3S le e 6 ), 

Dq = cosh(2zg + f3S se e 6 ) + 46cosh(2;| - /3<5i 6 e 6 ) + d + aa 6 + a/a 6 , 



z« = - hi 



1 1 + 7?W Z (6) 

1 -7 7 ( 1 )*(6) 



+ ^ c (6)7 ? ( 1 ) z (6)-^ 6 e 6 + %^, (2.10) 



a = cxp(— (3s), b = cxp(— f3w), d = exp(—(3wi), a 6 = exp(— pS^Se), 



also 



71 f W_ £><*(6)' 



(a = x,y), 



where 



m? (6)= sinh(^i^ + /?<5 s6 e 6 )+dsinh &±± 

3 Z Z 

. , A? (6) , a . , AJ(6) 
iflfle sinh 3V ; ±— sinh 4V y + 
2 a 6 2 



±smh( 5 V ; 



v4 Q (6) 

/3<5i 6 e 6 ) T sinh( 6 v - i35 w e 6 )- 



+ sinh(^^ - (35 16 s 6 ) + sinh(^^ - /W 16 e 6 ) 

m?(6)= sinh(^^ + /?<5 s6 £ 6 )-dsinh ± 
4 2 2 

±aa 6 sinh 3 v ; =F~ sinh 4 v y + 
2 ag 2 



+ b 



sinh(^^ - (3S 16 e 6 ) + sinh(^^ - (3S 16 e 6 )± 
± sinh(^^ - /W 16 e 6 ) T sinh(^^ _ p Sl6 s 6 ) 



£> Q (6)=cosh( 



Af(6) 



-<isinh 



A£(6) 



, A? (6) a , A2(6) 

+aa 6 cosh ^ y H cosh 4 „ ' + 

2 a 6 2 



cosh(^-^ - /3(5 16 e 6 ) + cosh(-^p 



(2.11) 



(2.12) 



- pS 16 e 6 )+ 



+ cosh(^^ - /W 16 e 6 ) + cosh(^^ - (38 16 e 6 ) 
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Here the following notations are used 



and 



'16 



"26 



"46 



'16 



"36 



"26 



A?(6) = z?(6)±z2{6) + 4(6) ± 4(6); 
A? (6) = 4(6) ± 4(6) - 4(6) T 4(6); 
Af (6) = ±4(6) + 4(6) T 4(6) + 4(6); 
(6) = 4(6) ± 4(6) + 4(6) T 4(6); 



1, l + r?P(6) 
^ 1 - jji (6) 



+ /Ji^(6)+/J«^ x (6 

1 1 + 7^(6) a ( 

z i - m ( 6 ) 

+ /jrf*(6)+/?rf x (6 



1, l + ^2 1)a: (6) fl 

+ /3^^ 1)a; (6)+/3^r7i 1)a; (6 

1 1 + ^(6) , . 
z 1 - ?74 (6) 

+ /3^ 3 ^ 1)a; (6)+^ 1 r7i 1)a; (6 



1 1 + ^(6) 

2 i _„(!)»/ 



■ /^lW 

1-^(6) 
+ /3^ 2 ^ 1)! '(6)+^ 2 r 7 | 1)! '(6 

1, l + ^ 1)y (6) , 
2 1-^(6) 

1 1 + ^(6) 



/3^ ± cos 7 Si 
f3i>e£6 + 2 ' 

1)x (6)+/3^^ 1)a: (6) + 
cos^E- 



■ Pi>eS6 

1)x (6) + /3^^ 1)a: (6) + 

/3^ ± sin7 J Bi 
We£e 2 : 

1)x (6)+^ 2 ^ 1)a: (6) + 

/3^ ± sin7 J Bi 
- /3^6£6 + ^ : 

1)y (6)+^ 1)?/ (6) + 



i (2.13) 



(6)+/?^ 1)!/ (6) 



sin7i^ 2 ,„ „ 
-/fr/> 6 e 6 + / ^ ± o r ,(2.14) 



1 )«(6)+/?^ 1)i '(6) + 

/3^ ± cos7£ 2 
- /3Ve£6 » , 
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\ In 1 + + /W^W + PwgtoQ) + 

z 1 — 774 y (6) 



'46 o *" . (l)y 



3 Elastic, piezoelectric, and dielectric 
properties of the KD2PO4 type crystals 
under the mechanic stress erg 

Influence of the stress ae on the characteristics of the KD2PO4 type crys- 
tals will be considered using the thermodynamic potential per primitive 
cell, which within the four-particle cluster approximation reads 

9ie(6) = |4 e| - ve° 36 e 6 E 3 - V - X f 3 Ej + 2Tln2 + (3.1) 
+ 2i/ c [t/ 1 > z (6)] 2 - 2Tln[l - (?7 (1) *(6)) 2 ] - ZTlnDg - vc7 6 e 6 . 

Let us mention the equation for the extremum of the the rmo dynamic 
potential g\E with respect to rf 1 ^ coincides with Equation 
From the thermodynamic equilibrium conditions 

v V 9e 6 J E3 r76 v V dE 3 J ae 

we obtain 

66 6 36 J -r _ VDI vD% vD% ' 

^ = e° 6£6 + x^3 + 2^^, (3.2) 
where 

M a6 = aa 6 -— , M s6 = sinh(2z z + /3<5 s6 £ 6 ), M i6 = 46sinh(zg -/3<5i 6 e 6 ). 
a 6 



From (3.2) we find the electric field 

E 3= _ h o 6£6+kl o(p 3 _ 2 !i!!my (3.3) 

wher e h 3 Fl = e 36 /x%3, k 33 = l/x§3- Substituting the expression (|3.3[) 
into ( |3.2|) , we get 

„ fcO f P 9 ^™ z (6^ , 4^ 6 m z (6) 

°6 - c 66 e 6 - h 36 f 3 - 2 — 



v v D% ) v D% 

-^M a6 - ^M s6 + ^M 16 , (3.4) 
vDl vDl vDl bl v ' 
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where c™ = eg + e° 3e h° m - 



Expressions ( p. 3] ) and (3.4) can also be obtained from the conditions 
df\ 1 / df\ 

whereas the free energy / is 

f(6)=g2E(6)~vP 3 E 3 . (3.5) 



Substituting the expressions ( |3.l|) and ( |3.3|) into fl3.5|) , we obtain 

/(6) = J^ 2 . - vh° m e 6 P 3 + - £*& ( 2 |^(6)) 2 + (3.6) 

+ 2Tln2 + 2j/ c [t7W 2 (6)] 2 - 2Tln[l - (?y (1)2: (6)) 2 ] - 2Tln£>f. 

Let us study now the influence of the stress a§ on piezoelectric, di- 
electric and elastic properties of KD2PO4. 

From the expression for the mean value of quasispin (<J q f) (2.9) it 
follows that 



de 6 ) Ea Dq- 2(f2^e ' 
where 

x e = cosh(2z 6 + (3S s6 s 6 ) + 6cosh(z 6 - (35i 6 e 6 ) - ?7 (1) (6)ra(6), 
r 6 = 5 s6 cosh(2z 6 + f3S se e 6 ) - 2b5i 6 cosh(z 6 - f3Si 6 e e ) + 
+j/ 1 >*(6) [-5 s6 M s6 + 6 a6 M a6 + 5 U M 16 ] , 

= 1 - (t?( 1 )(6)) 2 + PVc ' 06 = ~ 2 ^ 6 + r& - 
Hence, the coefficient of piezoelectric stress e 3 § is 

da 6 \ fdP 3 \ _ , 2fi 3 f39 6 



In the paraelectric phase at o§ — 0, the coefficient of piezoelectric stress 
e 3 § equals 



. _ 1 M3 2/3[-2(l + b)ip 6 + 6 s6 - 2b6 16 ] 
336 " 636 + v -l + 2b + 2a + d-2/3v c (l + by 



Differentiating the expression for polarization ( |3.2|) with respect to 
the field E 3 at constant strain Sq we obtain dielectric susceptibility of a 
clamped crystal 

V ^3 / e , v T D 6 - 2>c%tp e ta 6 
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Using the relations 



de 6 

n 

da 6 



= /}0 6 + /}/i 3 X 6 n — 
P 3 V ^6 



P 3 



^3 



9cr 6 



6/ _B 3 



where = 



9cr fi 



is the crystal compliance at constant electric field, 



we obtain the constant of piezoelectric stress h 3 g: 



h 36 = - 



fdE 3 



_( 9oq_ \ _ e 36 



6/p 3 



8ft/, 



X33 



and the coefficient of the piezoelectric strain 



^36 = e 36 s 



66- 



(3.10) 



(3.11) 



From (3.2) and (3.4) we get the expressions for the elastic constant 



at constant field 



' 0G 



4^r 6 2 



■60 



El 



u L> 6 - 2x 6 <^g vTD e (D 6 - 2x 6 <^) 



5 2 6 cosh(2z 6 +/35 s6 £ 6 ) + <^ 6 (ac 6 H — ) + 4M 2 6 cosh(z 6 -/3(5i 6 e 6 ) 

c 6 



2/3 2 
+ — p- [-<J S 6-M^s6 + S ae M a 6 + Si 6 M w ] ; 
vD 6 

and for the elastic constant c£ at constant polarization 



-66 — °66 



e 36 /i 



36- 



(3.12) 



(3.13) 



In the paraelectric phase at erg = when 7/^(6) = 0, the renormalizcd 
elastic constant for an unstrained crystal is 



L 66 



„E0 



4^6 



-2(1 + b)^ 6 + 5 s6 



vT vT[-l + 26 + 2a - 2/3z/ c (l + b)] 



(3.14) 



4^6 2 

vD+T vT(l + 4b + 2a) 



(<5 2 6 + 2ad 2 a6 + S 16 ). 



Hence, the relations ( p.2j) — (3.4) take the form 

°6 = cf 6 e 6 - e 36 £ 3 , a 6 = c£ 6 e 6 ~ h 36 P 3 , 
P3 = e 3 6£6 + X 33 E 3 , E 3 = -/i 36 e 6 + fcf 3 F 3 , 



(3.15) 
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where /c| 3 = 1 /X33 is the inverse dielectric susceptibility. These relations 
can be also obtained within the phcnomcnological approach, but we 
preferred to have the microscopic expressions for cf e , e 3 §, x 33 t h. 



36- 



From the system of equations (3.15) at E 3 = const 



■ 66 



'36 



6/ E 3 

de 6 
da e 



h>36 

+ k 



dP 3 
da 6 

dP 3 



= 1, 



E 3 



6/ B 3 



<9(7 6 



= 



6 / E 3 



we find the expressions for the coefficient of piezoelectric strain 

_ ^36 636 

£3 



^36 = 



dP 3 
da 6 



„P ue 



h 2 

"36 



-66 



36 



and compliances at constant field 



A 66 



33 



da 6 J 



„P he 

<-m K, 33 



h 2 
"36 



-66 



e 3 (ih 



36 



1 

' 66 



From (3.15) at P 3 = const we get 
de 6 



L 66 



e36 



6/P3 



<9er 6 / 



e36 



X33 



dE 3 \ 
da 6 J 

dE 3 \ 
da 6 J 



p 3 



= 0; 



P 3 



hence the constant of piezoelectric strain is 



536 



dPh 
daa 



e 3 e, 



'36 



P3 



X33 



-66 



e 3 &h 



66 



'36 



' 66 



and compliance at constant polarization is 



s 66 



_ [dee \ 1_ 



\9a 6 J 



-66 



66 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



Differentiating the system ( 3.15| ) with respect to the field E 3 at con- 
stant stress, we get 



(da 



-66 



\dE 3 



em = 0, 



-e 3 6 



Since 



'36 



\0E 3 

( d £ 6 
\dE 3 



fdP 3 



\dE 3 



= X33- (3-21) 
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then the dielectric susceptibility at <tq = const is 



[^jf] =\!u <:«<1:«, (3-22) 



E 



Hence, we obtained the microscopic expressions for Cg 6 , e36, and X33; the 
other characteristics are expressed via them. 

Let us calculate now the transverse dielectric susceptibility of the 
KD2PO4 type crystals in the presence of stress a§ defined as 

X11 = 7. 



v \ dE l dEi J 

/ij_sin7 / dri2^ x {&) dr)^ x (6) 
v I Mi + dEl 



(3.23) 



Using the relations ( [2.11 )-(2.14) we obtain the following system of equa- 
tions 

/3/z_L (cos 7^ - M q6 sin 7) = 

- & - m - M "*- m 

f3fi ± (M a6 cos 7 - sin 7x5) = 

d U^(6) - 4^(6)) d U^(6) - ^(6)) 
= (D 6 *^ 6 ) — M rf6 — , 

where the following notations are used 
a 



Xq =aa 6 H h2fecosh(z 6 -/3(Si 6 e 6 ), tp% = m +/3(i/ 1 -v 3 ). 

a 6 1 — iv (6)) 

In the result 



xh 



T (D 6 - <^<) 2 - (<^M a6 ) 2 
- ^4-7 ^ sin 2 7 . (3.24) 

4 Influence of stress cr 6 on thermal proper- 
ties of the KD2PO4 type crystals 

Molar entropy of the deuteron subsystem of KD2PO4 type crystals under 
stress o§ reads 

S 6 = - R f^ 



V dT 



P 3 ,£6 



i?|21n2 + 21n[l - (77^^ (6)) 2 ] + 2hiZ) 6 
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4T^ 7 «(6) + ^ ) (4.1) 



D 

where R is the gas constant, and 

^ = -^(^ (1) (6)-V 6 £ 6 ), (4.2) 

,,w , ,wi ( a \ e ( a \ 5qEq 
Mq — 46— cosh z 6 + a— + ac§ H — + ac 6 



T T \ cqJ T \ cqJ T 

The molar specific heat of the deuteron subsystem of KD2PO4 type 
crystals can be calculated by differentiating the entropy 

AQ = t(^J =AQ+q 6 p a 6 , (4.3) 

where AC| is the molar specific heat at constant strain 

ACI=q^+qlpl (4.4) 



Using the relations (4.1) we obtain 



= ^ |2T^[2x 6 T^ + 2(q 6 - r,«(6)M 6 )] + N 6 - M| , 
98 - fit) = -^^{2x 6 T^ + [ ft - ^)(6)Af«]} 



-6. 



is the polarization heat at given 

% P = (||) p t = ^ {2T^(-2x 6 ^6 + r 6 ) - 2[q 6 - 7jW(6)M 6 ]^ 6 
M 6 

— Ae + -f^-(S s6 M s6 + <S o6 M o6 + S 16 Mi 6 ) 

is the deformation heat at given Pg, where 
1 



(e + (5 a 6£6) aa e + (e - S a6 e 6 ) h 46w cosh(zg - /3(5i 6 e 6 )+ 

»6 



+ wi 2 d - (6 s6 e 6 ) cosh(2z 6 + /3S as e 6 ) + 

+ (^i6£6) 2 4&cosh(z 6 - (38 16 e 6 ) + (5 16 e 6 w) 86cosh(2: 6 - /3S 16 e 6 ) 
qci = ^[-5 s6 e 6 cosh(2z 6 + f35 as s 6 ) + <5i 6 e 6 26cosh(z 6 - /3<5i 6 £ 6 )4 
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+ 2bw cosh(z6 — /3<5i6£6)] , 
A 6 = Y [S%£e sinh(2z 6 + j35 as e & ) + Sj 6 e e 4b cosh(z 6 - /3<5i 6 e 6 )+ 

+S a6 [aa 6 (e + S a6 e 6 ) - — (e - <5 Q6 e 6 ) ] + (5i 6 we 6 46sinh(z 6 - f35 16 e 6 ) 



In (4.3) and (4.4) pg = (dPs/dT)^ Ea is the pyroelectric coefficient, and 
c*6 = (dee/dT) a is the thermal expansion coefficient. 
From relations ( 3.1 5| ) we get 

K=P6 + e 36ae, (4.6) 

where 

£ _ Ms 2 2^ 6 r^ + [g 6 -7/ 1 )(6)M 6 ] 

and the thermal expansion coefficient reads 

a e = ^ , (4.8) 

c 66 

where ^6 = {jjjr) e — Qe ^ s the thermal pressure. 



5 Discussion 



Before going into the discussion of the proposed in previous sections 
theory, let us note that, strictly speaking, this theory can be valid for 
a completely deuterated KD2PO4 crystal only, whereas the vast ma- 
jority of experimental data concerns the crystals, deuterated partially 
(see [l2| [l3[ |l4], |fl| H). Nevertheless, as the relaxational character 
of dielectric dispersion in K(Hi_ x D :Z ;)2P04 crystals implies, tunneling 
effects in them are most likely suppressed by the short-range correla- 
tions between hydrogens 0, |2(], So at least in the case of high 
deuteration, we are allowed to neglect tunneling and apply the theory 
to partially deuterated crystals as well. Hereafter we shall consider the 
crystal with nominal deuteration x = 0.89 with the transition tempera- 
ture Tco = 210.7 K, for which experimental data for the relevant elastic, 
piezoelectric, and dielectric characteristics are available. 

As the values of the theory parameters e and w we choose those found 
in |2^, |2^| , which at <jq = provide a satisfactory agreement of the cal- 
culated curves with the experimental data for spontaneous polarization, 
specific heat, static and dynamic dielectric permittivities of the crystal. 

To determine the deformation potentials and "seed" quantities, we 
use the experimental data of Jl5|, [l6) for the temperature dependences of 
the coefficient of piezoelectric strain d,36, dielectric susceptibility X33, and 



16 



compliance sf 6 at as — 0. Using the known relations between dielectric, 
elastic and piezoelectric characteristics and having the values of X33J 
and Sq 6 we can calculate the piezoelectric constants e^e, h^e, <?36, elastic 
characteristics qf 6 , cf 6 , Sq 6 , and dielectric susceptibility x§3- 

Analogous calculations were carried out also for x = using the data 



of |l2j, O, 14 . Let us mentions that the experimental values of d^e or 
Cg 6 do not agree with the results presented in PH . 

The deformation potentials i/jq and ^s6 , and the parameter of the long- 
range interaction v were found from the condition that the transition 
temperature at = was Too = 210.7 K, and that the best description 
of the e36, s|g, and £33 temperature curves was obtained. Let us note 
that we get the best fit if the following constraint 

w 

S s6 = 325 + 2^ 6 exp(- — ) (5.9) 
J- CO 

is held. By changing tpe ( or ^se) we can slightly alter the theoretical 
slopes dh^/dT and dg^^/dT at T > Tco- So, fitting the theoretical 
temperature dependences of the piezomodules h^e and 336 to experi- 
mental points, we choose the values of ipe and S s e. It should, however, 
be noted that we can vary (5 s6 rather strongly (from up to the adopted 



in this paper value), provided Eq. 5.9 is obeyed, and still obtain a fair 
description of /136 and g^Q. 

The adopted value of the long-range interaction parameter is much 
lower than the one used in the theories where spontaneous strain is not 
taken into account [^2[ Low v allows us to describe the para- 

electric Curie constant and ferroelectric spontaneous polarization with 
the same value of the effective dipole moment ^3 describes response of 
heavy ions to deuteron ordering) , strictly determined by the experimen- 
tal saturation polarization. 

The values of 5 a 6 and were chosen by fitting to experimental data 
the calculated sf 6 (T) and c§ 6 (T) dependences. Moderate values of <5i6 
(splitting of the single-ionized deuteron configurations) do not percep- 
tibly affect the calculated curves, and its larger values only make the 
fitting worse. Hence, for the sake of simplicity we neglect this splitting, 
taking Sie = 0. 

The "seed" e^jg and X33 are merely the high temperature limits of 
experimental temperature dependences e^e and x| 3 . 

The adopted values of the theory parameters are presented in Table [IJ. 

To find the order parameter rp~> and strain Eq we minimize the ther- 
modynamic potential <7ib(6) with respect to 77 W and determine e§ from 



Eq. fly) 



Let us demonstrate how the presented theory desribes the physical 
characteristics of the crystal related to the strain £g and polarization 
P3 at erg = 0. As one can see, the theoretical results are in a good 
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Table 1: The theory parameters. 



£ W V c 


tf>6 S 3 6 8 a e, Sis 


C i,°-10- 10 &/v eg 6 xS J 


(K) 


(K) 


(dyn/cm^) (/iC/cm^) (esu/cm^) 


88.3 778 35.976 


-250 1750 -187.5 


7.0 6.21 0.42 ■ 10 4 0.4 



quantitative agreement with the experimental data in the paraelectric 
phase. 

In fig. H we plot the dependence of the compliance sf 6 and elastic 
constants cf 6 and Cq 6 of the K(H .nD .89)2PO4 crystal on AT = T — Tc 
at erg = 0. Here we also presented the experimental data for x = 0.89 
|U and x = 0.00 Jn|,[l4j. As 

one can see, the theoretical results are in 
a good quantitative agreement with the experimental data of |jl6| in the 
paraelectric phase. The compliance sf 6 at T — > Tq has an anomalous in- 
crease, whereas the elastic constant Cq 6 vanishes at the Curie point. The 
elastic constant Cg 6 at T < Tq is nearly constant with temperature, only 
has an about 6% decrease at T = Tq , and slightly increases in the para- 
electric phase. That accords with the conclusions of p5| , where it was 
shown that the ferroelectric phase transition in KH2PO4 did not affect 
the magnitude of Sq 6 , and hence of Cq S . According to the data of | Q , 
at x = 0.0 the elastic constant cjf 6 slightly decreases with temperature 
in the paraelectric phase. 




o 
o 



I , 1 , 1 , 1 , 1 o I , Q , 1 , 1 , 1 

-25 25 50 A7XK) "25 25 50 AT(K) 



Figure 2: Temperature dependence of compliance sf 6 : • - |lq|, A - sf 6 = 
I/cm H| and elastic constants Cg 6 an d cf fi : A - c<f 6 = l/s 66 [jl6), c<f 6 = 
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Temperature dependence of the coefficient of piezoelectric strain 
and coefficient of piezoelectric stress e^e of the K(H .iiD .89)2PO4 crys- 
tal at (76 = is shown in fig. || A satisfactory agreement of the presented 
theory with the experimental data of 15 for and with the data ob- 
tained from the formula e^e = <^36 HTsfeQ is observed. At T -> T c 
coefficients d^e and e^e sharply increase. In the ferroelectric phase, the 
calculated coefficients d^e and e^e sharply decrease with a rate much 
higher than in the paraelectric phase. 




Figure 3: Temperature dependence of coefficient of piezoelectric strain dm 
(left, experimentaf points taken from |p[ - • (x = 0.89), 0, - o (x = 0.00), 
|lf| ] - □ Ix. = 0.00)) and coefficient of piezoelectric stress (right, A - 




e:i6 = 



Temperature dependence of the calculated constant of piezoelectric 
stress /136 and constant of piezoelectric strain 1736 at <tq = is shown 
in fig. [I] along with the obtained using the experimental data of |l4], 
|l5[ values of h 36 and 336 for x = 0.89 and the data of |llj] for x = 
0.00. The piezoelectric constants h 3e or g 36 do not have singularities at 
the ferroelectric transition; therefore they are called "true" piezoelectric 
constants of the crystals. 

In fig. U we plot the temperature dependences of the dielectric per- 
mittivities of a free (efg) and clamped (e| 3 ) crystals at er 6 = along 
with the experimental data of Jl5| for 633. The "experimental" data 
for the dielectric permittivity £| 3 were calculated using the relation 
e 33 = £ 33 ~ 47re 36 d3 6 . 
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Figure 4: Temperature dependence of constant of piezoelectric stress hsa (left, 
A - /136 — d^mM/ Ss 6 X33WMi Q ~ vM) an< i constant of piezoelectric strain #36 
(right, A - fl 3 6 = dse/xh if, □ ~ @)- 

The presented graphs indicate that one can describe the tempera- 
ture peculiarities of these dielectric, piezoelectric, and elastic properties 
of KD2PO4 attributing those pecularities to deuteron subsystem only, 
with the heavy ions lattice counterpart considered as background and 
temperature independent one. However, since we cannot unambigously 
set the values of ipe and S s e, we are not in position to determine what 
are the relative weights in these peculiarities of piezoelectric coupling, 
described by the parameter ip 6 , and of the short-range up-down deuteron 
configurations splitting, induced by strain o§ and described by the pa- 
rameter (5 S 6 . 

In fig. ^ the dependences of the thermodynamic potential g\E on the 
order parameter t/ 1 ' at different values of temperature and stress o§ are 
presented. The behavior of g\E at zero stress ag is usual for the first 
order phase transition: a little below the transition point g±E has three 
minina - one at t/ 1 ) = and two symmetric ones ir]^ 0. The last two 
are of the same depth and lower than the one at rf 1 ^ — 0. At T = Tq 
all minima are of the same depth (criterion of the phase transition), 
and at T > Tc the central minimum becomes the deepest. At higher 
temperatures the minima of g\E at 77W ^ disappear. 

Under stress the dependence of giE(v^) becomes asymmetric; 
the lower values of the thermodynamic potential are at those values of 
the order parameter, the sign of which coincide with the sign of the 
stress. 

In fig. we plot the dependences of the solutions of the equation 
for the thermodynamic potential g\E extremum and the corresponding 
values of g\E on temperature at different values of stress ag i n the vicinity 
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Figure 5: Temperature dependence of dielectric permittivity of a free crystal 
£33 (left, ■ - pa]) and dielectric permittivity of a clamped crystal £33 (right, 

A-£l 3 = £ 3 CT 3 -^ 

°66 



of the transition point. 

Equation (2.9) may have up to five different solutions, three of which 
correspond to minimuma, and two correspond to maxima of the ther- 
modynamic potential. One of the minimuma - the central one - is at 
small values of the order parameter, and the sign of rp~> coincides with 
the sign of the stress <tq. The two other minima, if they exist, are nearly 
symmetric (absolutely symmetric at ctq — 0), and, as we have already 
mentioned, that minimum, the sign of which coincides with the sign of 
the stress is deeper. 



Let 77^ < 0, rj^ 



0, and ?7g 1 ' 1 > be the possible solutions of 



9lE 



9 1E 




Figure 6: Dependences of the thermodynamic potential g\E on the order 
parameter rf 1 ' at different values of stress erg. 
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Figure 7: Temperature dependences of the thermodynamic potential and 
solutions of equation for the order parameter at different values of stress u 6 : 
a) bar, b) 40 bar, c) 77 bar, d) 100 bar. o and • correspond to maxima and 
minima of the thermodynamic potential, respectively. 
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Equation ( |2.9| ), corresponding to the minimuma of g\E- As one can see 
(fig. |), at a 6 > 

9ie(v { 3 1} ) <9ie(t^) atT < T c , 

9ie{<1z^) = 9ie{V2^) at T = Tc, (the phase transition criterion) 
9ie(4 1) )>9ie(4 1) ) at T>T c . 

The temperature of the first order phase transition at which the 
branches of the thermodynamic potential intersect, increases with stress 
tJ6, and the values of r)^ at the transition point increase, whereas 773 
decrease. A decrease in the magnitude of jump in the order parameter 
( 5 7 y( 1 ) — _ yy' 1 ' means that the order of the phase transition is moving 
towards to the second order. At certain stress <7g Srj^ turns to zero - 
there is a critical point where the second order phase transition takes 
place (at T*). Further increase in stress smears off the phase transi- 
tion and results in continuous and smooth temperature dependence of 
the order parameter. Such behavior is typical for the first order phase 
transitions in ferroelectrics to which the electric field conjugate to spon- 
taneous polarization is applied p6fl . 

Corresponding Tc — ere phase diagram is depicted in fig.^[ Only 
the stable phases corresponding to absolute minima of thermodynamic 
potential are showed. The following temperatures are indicated: Too — 
210.7 K is the temperature of the first order phase transition at o§ = 0; 
T* = 211.6 K is the temperature of the second order phase transition 
at the critical point at stress a% = 75 bar; To = 210.8 K is the Curie 
temperature of a free crystal ag = 0, at which the compliance Sq 6 and 
the longitudinal static dielectric permittivity £33 diverge. As one can 
see, the diagram is symmetric with respect to a change — * —&6- An 
increase in the transition temperature with stress erg is practically linear 
with the slope dTc/d\ae\ = 21.6 K/kbar, which is much faster with than 
with hydrostatic or uniaxial p = —03 pressure [g7| . 

The observed "Y" -shaped form of the phase diagram is not unique 
but typical for the systems in fields conjugate to the order parameter. 
The Tc — E3 phase diagram of KH2PO4 of the same topology was ob- 
tained within the phenomenological approach without invoking any mi- 
croscopic model by Schmidt p8| ; an increase in the transition temper- 
ature with £"3 and cease of the transition at certain critical EZ were 
observed experimentally in KD2PO4 by Sidnenko and Gladkii p9fl . 

Let us note that in order to observe the dependence of temperature 
phase transition on stress o§ experimentally, one should apply an ex- 
ternal stress to the paraelectric crystal and cool it below the transition 
point. In absence of external fields or stresses, crystal at the transition 
can with 50:50 probability have positive of negative spontaneous strain. 
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Figure 8: Tc — ere phase diagrams of a KD2PO4 crystal. 



However, since in the presence of the stress o§ the minimum of the ther- 
modynamic potential at 7/ 1 ) of the same sign as that of a§ is deeper 
than the opposite minimum, the sign of the strain below the transition 
point will coincide with the sign of stress (similarly, as the direction of 
polarization in a crystal coincides with the direction of external field). 
The same picture will take place if the stress is applied to spontaneously 
polarized and strained crystal, provided that this stress induces the strain 
of the same sign as that of the spontaneous strain. However, if stress 
is applied to spontaneously polarized and strained crystal such that the 
signs of the stress induced and spontaneous strains are opposite, the sys- 
tem appears in a local minimum of the potential (for instance, at o~q > 

9ie(t)i) > 9iE(Wa))) that is, in a metastable state. To predict theoret- 
ically when the transition to a stable state (which requires switching of 
polarization and strain) is impossible. Temperature of such a transition 
will essentially depends on the experimental conditions. Nevertheless, we 
can maintain that this temperature is not higher than the temperature 
of absolute instability (the point where the local minimum disappears), 
and, as one can see at Figure [?], the temperature of absolute instability 
decreases with stress cr 6 . 

Let us consider now the dependences of the values of the order pa- 
rameter corresponding to minima of the thermodynamic potential on 
stress erg at different values of temperature in the vicinity of the tran- 
sition point, depicted in fig. 0. At T < Too, the branches of the ther- 
modynamic potential giE(Vi < 0) an d 9ie(ti^ > 0) intersect, and 
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the solutions of Equation ( |2.9| ) rj[ ' and 773 can exist simultaneously. 
Experimentally, on changing stress ag a regular hysteresis loop P3 — ag 
should be observed. 

At Tqo < T < T* , only one of the non-zero minima and the central 
minimum of the thermodynamic potential (tj^ and r/^ or ^3 an d 
772 ) can coexist. Experiment should reveal a double hysteresis loop. At 
temperatures higher than the critical T*, the dependences <?i_e((76) and 
^H ^) are smooth, and no jump in the order parameter is observed. 
Let us note that such a sequence of the hysteresis loops - a single loop 
at T < Tqo, a double loop at Tqo < T < T*, and a gradual change 
at T > T* - was obtained yet by Merz for BaTi03 in electric field fl30f 
and by Sidnenko and Gladkii pj| for KD2PO4 in the field E3, and is, 
apparently, typical for the phase order phase transitions in ferroelectrics 
under fields conjugate to the order parameter. 

In fig. [l(] we present the temperature and stress dependences of strain 
£6 at different values of stress eg and different values of temperature. 
Fig. [n] illustrates the analogous dependences of polarization P3. As one 
can sec, these quantities exhibit the identical variation with temperature 
or stress. The spontaneous strain £g and polarization P3 at <7g = in 
the ferroelectric phase slightly decrease with temperature and jump to 
zero at T = Tq ■ When stress a§ is applied, the magnitudes of strains £6 
and polarization P3 increase, whereas the jumps A£6 and AP3 decrease; 
at (76 = erg and T = Tq, the jumps vanish. At erg 7^ £6 and P3 
differ from zero above the phase transition and slightly decrease with 
temperature. At Tq < T < values of £6 and P3 increase with stress 
<76) having upward jumps at stresses corresponding to the curve of phase 
equilibrium with A£6 and AP3 decreasing with stress and vanishing at 

^6 = O-g*. 
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Figure 9: Stress dependences of the thermodynamic potential and solutions of 
equation for the order parameter at different values of temperature: a) 209 K, 
b) 210.5 K, c) 211 K, d) 212 K. T co = 210.7 K. o and • correspond to maxima 
and minima of the thermodynamic potential, respectively. 
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Figure 10: Dependences of strain eq on temperature at different values of 
stress a 6 (bar) (left): 1 - 0; 2 - 40; 3 - 40; 4 - 100; 5 - 150 and on stress a 6 
at different values of temperature T(K) (right): 1 - 209.0; 2 - Too = 210.7; 3 
- 211; 4 - 211.3; 5 - T* = 211.57; 6 - 211.8. 
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Figure 11: Dependences of polarization P3 on temperature at different values 
of stress a 6 (bar) (left): 1 - 0; 2 - 40; 3 - a% = 75; 4 - 100; 5 - 150 and 
on stress <jq at different values of temperature T(K) (right): 1 - 209.0; 2 - 
Tco = 210.7; 3 - 211; 4 - 211.3; 5 - T* = 211.57; 6 - 211.8. 

Further increase in stress slightly raises up the strain e 6 and polar- 
ization P 3 . At temperatures T <Tq and T — Tq > 3 K values of e 6 and 
P3 exhibit a nearly linear increase with stress ag. At changing the sign 
of the stress o§ , the absolute values of the strain £ 6 and polarization P 3 
do not change, but become negative. 

Let us discuss the theoretical results for the stress influence on the 
physical characteristics of the K(H .iiD .89)2PO4 crystal. The ma- 
jor stress eg effects on the physical characteristics are related to the 
changes induced by this stress in the character of the phase transition 
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and, therefore, essential only in the vicinity of the transition point; at 
\T — Tc\ > 5 K, these effects are negligibly small. 

In fig. |l2] we plot the temperature and stress dependences of the 
elastic constant Cgg at different values of stress a§ and temperature. As 
one can see, the elastic constant qf 6 change with temperature and stress 
exactly as the strain £6 and polarization P3 do. 

c 66 p , 10"dyn/cm 2 c 6 /, 10 1 'dyn/cm 2 



6.75 - 




210 21 1 212 213 T(K) 50 100 150 &(,(ba.l) 



Figure 12: Dependences of elastic constant c^ 6 on temperature at different 
values of stress cr 6 (bar) (left): 1 - 0; 2 - 40; 3 - ol = 75; 4 - 100; 5 - 150 
and on stress <T6 at different values of temperature T(K) (right): 1 - 209.0; 2 
- Too = 210.7; 3 - 211; 4 - 211.3; 5 - T* = 211.57; 6 - 211.8. 

Temperature dependences of compliance sf 6 (fig. |l3|), coefficient of 
the piezoelectric strain e?3g and stress e^Q (fig- |l|), an d of the dielec- 
tric permittivity e| 3 (fig. |l^) at different values of stress og are similar. 
The maxima of these characteristics shift to higher temperatures with 
stress. At er 6 < <7g, these characteristics reaches their maximal values 
at the transition temperatures. At eg = Og T = T£ the second order 
phase transition takes place accompanied by a sharp increase in these 
quantities. 

The temperature curves of the constants of piezoelectric stress 536 
and piezoelectric strain g^Q (fig. |l6| ) at different values of stress oq are 
similar to each other. At ae < Cg, minor jumps in h^e and f/36 are 
observed, which vanish at erg = erg. At erg > erg /i 3 g g^Q increase with 
temperature, tending asymptotically to the values of h^e, g^e at o§ = 0. 
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Figure 13: Temperature dependence of compliance Sqq at different values of 
a 6 (bar): 1 - 0; 2 - 40; 3 - a* & = 75; 4 - 100; 5 - 150. 
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Figure 14: Temperature dependences of coefficient of piezoelectric strain 
and coefficient of piezoelectric stress es6 at different values of stress a% (bar): 
1 - 0; 2 - 40; 3 - a% = 75; 4 - 100; 5 - 150. 
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Figure 15: Dependences of dielectric permittivity £§ 3 on temperature at dif- 
ferent values of stress at, (bar): 1 - 0; 2 - 40; 3 - a% = 75; 4 - 100; 5 - 
150. 




Figure 16: Temperature dependence of constant of piezoelectric stress /136 
and constant of piezoelectric strain at different values of stress ae (bar): 1 
- 0; 2 - 40; 3 - erg = 75; 4 - 100; 5 - 150. 
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The stress dependences of the compliance sf 6 , coefficients of piezo- 
electric strain d^e, piezoelectric stress e36, dielectric permittivity e| 3 at 
different values of temperature are presented in figs, [l7|, [l^, 19. 



At temperatures below To, these characteristics decrease with stress 
(76, decreasing being the stronger, the closer temperature is to Tco- At 
Tea < T < Tq, the mentioned characteristics increase with stress ag, 
having downward jumps at stresses corresponding to the curve of phase 
equilibrium, with the jumps vanishing at T = and eg = a%- The 
further increase in stress lowers down sf 6 , d^e, e^e and e| 3 . At T > 
Tq, these quantities exhibit a gradual increase with stress, reach their 
maxima, and decrease on further increase in stress. As temperature 
rises, the maximal values of the quantities decrease and shift to higher 
stresses. At AT > 15 K, the values of s^, c^, e^a, and e| 3 are stress og 
independent. 



S 66 c , 10" n cm z /dyn 
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Figure 17: Dependences of compliance on stress a& at different values 
of temperature T (K): 1 - 209.0; 2 - T co = 210.7; 3 - 211; 4 - 211.3; 5 - 
T* = 211.57; 6 - 211.8. 
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Figure 18: Dependences of coefficient of piezoelectric strain d^e and coefficient 
of piezoelectric stress e36 on stress o§ at different temperatures T (K): 1 - 
209.0; 2 - T co = 210.7; 3 - 211; 4 - 211.3; 5 - T* = 211.57; 6 - 211.8. 




Figure 19: Dependences of dielectric permittivity e| 3 on stress aa at different 
values of temperature T (K): 1 - 210.0; 2 - T C0 = 210.7; 3 - 211; 4 - 211.3; 5 
-T* = 211.57; 6 - 211.8. 

In fig.^Ojwe depict the dependences of constants of piezoelectric stress 
/i36 and piezoelectric strain g^Q on stress ae at different values of tem- 
perature. At temperatures T < Tq the values of and g^e decrease 
with stress <j§. If Tq < T < Tq, the constants and g^ are slightly 
lowered down by stress. At stresses corresponding to the curve of phase 
equilibrium they decrease with a jump, and at a§ = <7g the jump mag- 
nitude turns to zero, and the further increase in stress gives rise to an 
almost linear decrease in and 536. At temperatures T > Tq and 
<736 exhibit a gradual decrease with stress ag 

Let us consider now the temperature and stress dependence of the 
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Figure 20: Dependences of constant of piezoelectric stress /136 and constant of 
piezoelectric strain on stress <T6 at different values of temperature T (K): 
1 - 209.0; 2 - Too = 210.7; 3 - 211; 4 - 211.3; 5 - T* = 211.57; 6 - 211.8. 



transverse dielectric permittivity en. Calculations were carried out at 
v x - f3 = 10K, fix = 2.79 • lCT 18 esu • cm, e^ = 10. In fig.|l| we depict 
the dependence of en on temperature at different values of stress erg and 
on stress at different values of temperature along with the experimental 
points for erg = 0. 

In the paraelectric phase, the experimental data are described by the 
theory rather well. As temperature rises, en gradually increases in the 
ferroelectric phase, has an upward jump at at the transition tempera- 
ture, the jump magnitude being equal to zero at erg = erg. At erg = 
the values of en slightly decreases with temperature in the paraelec- 
tric phase, whereas at erg ^ after the upward jumps it asymptotically 
tends to the values of en at erg — 0. At erg — erg en exhibits only smooth 
increase with temperature. 

At temperatures T < Tq, the magnitude of en is lowered down by 
stress erg. If Tq < T < Tq, en slightly decreases with stress, then at 
stresses corresponding to the curve of phase equilibrium it has a down- 
ward jump, the size of which Aen is zero at erg = Cg. At further increase 
in stress erg, en decreases linearly. At T > Tq the permittivity en ex- 
hibits a smooth decrease with stress erg. 
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Figure 21: Dependences of dielectric permittivity sfj on temperature at dif- 
ferent values of stress cr 6 (bar) (left): 1 - 0; 2 - 40; 3 - erg = 75; 4 - 100; 5 
- 150 and on stress erg at different values of temperature T(K) (right): 1 - 
209.0; 2 - Too = 210.7; 3 - 211; 4 - 211.3; 5 - T* = 211.57; 6 - 211.8. • - Q. 

In fig.^H we plot the dependence of specific heat of a deuteron subsys- 
tem of K(H .iiD .89)2PO4 on temperature at different values of stress erg 
and dependence of ACg on stress eg at different values of temperature. 



AC/, J/mol K 




r(K) 




50 a f (bar) 



Figure 22: Dependences of specific heat AC'q on temperature at different 
values of stress a 6 (bar) (left): 1 - 0; 2 - 40; 3 - a% = 75; 4 - 100; 5 - 150 
and on stress <T6 at different values of temperature T(K) (right): 1 - 209.0; 2 
- Tco = 210.7; 3 - 211; 4 - 211.3; 5 - T* = 211.57; 6 - 211.8. 

In the ferroelectric phase at erg — 0, the magnitude of ACg increases 
as temperature approaches Tq, has a downward jump at the transition 
point, and is nearly constant in the paraelectric phase. On increasing 
stress erg up to dg, the maximum of ACg increases, reaching its maximal 
value at erg = <7g, whereas the further increase in stress lowers down the 
maximal value of ACg and shifts it to higher temperatures. Under stress 
erg the jump ACg is smeared out. 
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At temperatures lower than Trj, the ACg increases with stress erg. 
At T < the magnitude of ACg first slightly increases with stress, 
and then at stresses corresponding to the curve of phase equilibrium it 
exhibits a sharp increase, diverging at T — Tq and erg = Cg ACg . The 
further increase in stress gradually decreases AC@ . At temperatures 
T > Tq, the specific heat ACg smoothly increases with stress, reaches a 
maximum and smoothly decreases. 

Temperature and stress dependences of the pyroelectric coefficient pg 
and of the coefficient of linear expansion ag are analogous to those of 
the specific heat ACg . 

6 Concluding remarks 

In this paper we present the microscopic model for a description of the 
stress ag on the phase transition, static dielectric, elastic, piezoelectric, 
and thermal properties of deuterated ferroelectrics of the KD2P04-type. 
Unlike hydrostatic or uniaxial p = —03 pressures, the stress ag lowers the 
symmetry of the high-temperature phase down to the symmetry of the 
low-temperature phase. As the electric field £3 , the stress erg induces in 
KD 2 P04-type crystals the strain e 6 and, due to the piezoelectric effect, 
polarization P 3 . 

In our previous papers we showed that an important role in depen- 
dences of the transition temperature and dielectric characteristics of hy- 
drogen bonded crystals of the KH 2 PC>4-family on pressures that do not 
change the system symmetry is played by the corresponding changes in 
hydrogen bond geometry, in particular, in the distance 5 between pos- 
sible deuteron sites on a bond. Most likely, the stress ag does not per- 
ceptibly affects S. Instead, it alters the angle between hydrogen bonds, 
perpendicular in an unstrained paraelectric crystal, and distorts the PO4 
tetrahedra, splitting thereby energies of deuteron configurations. An- 
other important mechanism of the stress erg influence on the phase tran- 
sition and physical properties of the KD 2 P04-type ferroelectrics is the 
piezoelectric coupling, giving rise to effective fields, action of which, for 
the symmetry reasons, is equivalent to action of an external electric field 
applied along the ferroelectric axis. 

As we show, the transition temperature increases with the stress erg, 
and the order of the phase transition tends to the second one. In the con- 
structed phase diagram, which has the same topology as the predicted in 
the phenomenologic approach Tq—E 3 diagram, there are two symmetric 
critical points where the second order phase transitions take place, and 
the curves of phase equilibrium terminate. The stresses above critical 
smear off the phase transition and smoothen the temperature depen- 
dences of polarization and strain. Correspondingly, an increase in the 
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stress eg raises up the peak values of the physical characteristics of a 
crystal having peculiarities in the transition point (the longitudinal di- 
electric permittivity, compliance sf 6 , piezomodules and e36, and the 
specific heat). These peak values are the highest at critical stresses, the 
higher stresses smoothen the temperature dependences of these charac- 
teristics and lower down the peaks. 

Taking into account the piezoelectric effect in the developed model 
allows one, at the proper choice of the theory parameters, to quanti- 
tatively describe the available experimental data for the temperature 
dependences of dielectric, elastic, piezoelectric, and thermal characteris- 
tics of unstrained KD2PO4. Further experimental studies are required 
to ascertain the values of the theory parameters and verify its predic- 
tions about the form of the Tq — <t& phase diagram and possible stress 
(76 effects on the physical characteristics of the crystals. 
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